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ABSTRACT 

We present the 2012 Hubble Ultra Deep Field campaign (UDF12), a large 128-orbit Cycle 19 HST 
program aimed at extending previous WFC3 /IR observations of the UDF by quadrupling the exposure 
time in the F105W filter, imaging in an additional F140W filter, and extending the F160W exposure 
time by 50%. The principal scientific goal of this project is to determine whether galaxies reionized 
the universe; our observations are designed to provide a robust determination of the star formation 
density at z > 8, improve measurements of the ultraviolet continuum slope at z ~ 7 — 8, facilitate the 
construction of new samples of z ~ 9 — 10 candidates, and enable the detection of sources uptoz~ 12. 
For this project we committed to combining these and other WFC3/IR imaging observations of the 
UDF area into a single homogeneous dataset, to provide the deepest near-infrared observations of 
the sky currently achievable. In this paper we present the observational overview of the project, 
motivated by its scientific goals, and describe the procedures used in reducing the data as well as the 
final products that are produced. We have used the most up up-to-date methods for calibrating and 
combining the images, in particular paying attention to correcting several instrumental effects. We 
release the full combined mosaics, comprising a single, unified set of mosaics of the UDF, providing 
the deepest near-infrared blank-field view of the universe obtained to date, reaching magnitudes as 
deep as AB ~ 30 in the near-infrared, and yielding a legacy dataset on this field of lasting scientific 
value to the community. 

Subject headings: Cosmology: observations — Galaxies: high-redshift — 



1. INTRODUCTION 

A fundamental quest of modern observational cosmol- 
ogy involves expanding the frontiers of knowledge about 
the formation of the first stars and galaxies at the ear- 
liest epochs of cosmic time, and determining their role 
in the reionization of the universe at redshifts above 7. 
This is also among the most challenging of observational 
regimes to explore, requiring depths up to 30th magni- 
tude (AB) or beyond, at count-rates that are thousands 
to millions of times fainter than the typical ground-based 
sky brightness per square arcsecond at optical and near- 
infrared wavelengths, respectively. 

The Hubble Space Telescope (HST) has played a unique 
role in these explorations of the very early universe, prob- 
ing these extreme depths by virtue of its combination of 
high angular resolution, and low sky background achiev- 
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able only from space. Significant investments of time 
in deep single-pointing surveys with HST have yielded 
a wealth of scientific results to date, from the original 
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Hubble Deep Fiel d dWilliams et al. 1996) and Hubble 



2000 



Deep Field South (Casertano et al.|l 



illiams et al. 



Lucas et al.][2003J, together with the 2004 Ultra 



Deep Field (UDF ^Beckwith et al.|2006| |Thompson et al. 
2005 1 , which has since become the centerpoint for deep 
follow-on im aging programs in 2005 ( PL: M. Stiavelli, 



described in Oesch et al. 2007 2009), as well as with 



the Wide Field Camera 3 lntrareTchannel (WF C3/IR) 



in 2009 (UDF0 9, PI. G. Illingworth, desc ribed in Oesch 



et al.||2010b|a 



and Bouwens et al. 



cintly in M 2 (UBFT2; Pi: R. El 



2011b I and mo st re 
described in Ellis 



et al.|2012 together with the present paper) Fj 

These surveys, in conjunction with wider, shallower 
HST surveys with the A dvanced Camera for Surveys 
ACS) including GOODS flGiavalisco et al.|2004| ), GEMS 
Rix et al.|[2004l, AEGIS (|Davis et al.][2U07]), C OSMOS 
Scoville et al ||2007| |KoeEio|ife^al.| |2007L WFC3 
Liarly Relea se science prog ram (Windhorst c t~al. | |20lT] ), 
CAND ELS (|Grogin et all 2011 1 [Koekemoer et aflpUTTt , 
BoRG ( |Trenti et_al.||2011|) HIP PIES ( |Yan et al.||2011b| , 
and CLASH |Postman ~et al.||2012 ), have transformed 
our understanding ot the early universe. There is now 
overall evidence for the mass build-up of early galaxies 
at z ~ 4 — 8 based on the evolution of the cosmic star- 



formation density (Giavalisco et al. 2004 Bouwens et al 



10 Note that there is a separate UDF imaging program, in the ul- 
traviolet with WFC3/UVIS (PL: H. Teplitz); this probes lower red- 
shifts complementary to the very distant universe discussed here. 
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Table 1 

UDF12 Observing Summary (HST Program ID 12498) 



Field 


Instrument / Camera 


Filter 


# Orbits 


# Exposures 


Exposure Time (s) 


UDF-main 


WFC3/IR 


F105W 


72 


144 


198,423 




WFC3/IR 


F140W 


.30 


60 


82,676 




WFC3/IR 


F160W 


26 


52 


71,652 


UDF-par2 


ACS/WFC 


F814W 


128 


256 


322,944 



20041 [20071 [Bunker et al.|200 4, McLu re et al 120061 [20091 
Yan et al12006[|2010[|Castellano et al.|2010[|Oesch et al'| 



2010bp . A wider variety of results have been obtained on 
the ultraviolet spectral slopes and stellar popu lations of 
these early star-forming galaxies at z ~ 7 — 8 ( Bouwens 
et~aL|[2009l [20101 [20121 |Q no et al.|[2010l |Bunker et al . 
2010[|Finkelstein et al.|2010[|2012b|aHYan et al.|2011b|a 
McLure et al.||2010i |2011[ [Grazian et al | |2011i |2012 



Bradley et al.||20121 |Dunlop et al.||2012a|b ) , particularly 



the deep 
which for the UDF 



concerning the extent to which these galaxies may or may 
not be able to account for reionization. Finally, tantaliz- 
ing discoveries of galax ies at z ~ 9 — 10 (Yan et al.| | 2010l 
Bouwens et~aLl|2011a[ |Zheng et al.||20l2p , z pTT^oe 
et al.|[20l^ ), and most recently up to z ~ 12 (|Fllis et al. 
2012p are only becoming possible by means o: 
est near-IR observations achievable, which fo: 
are described in this paper. 

The structure of this paper is as follows. The sur- 
vey overview and observational design are presented in 
£[2j followed by the description of the data processing 
and calibration in f|3j the presentation of the final data 
products in Sj4j and the summary in ^j5] Further details 
and current updates about the survey are provided at 
the project websitcp"| while all the final combined mosaic 
data products from the survey are bein g m ade publicly 
available as High-Level Science Product^] that are de- 
livered to the Space Telescope Science Institute archive. 
These images constitute a single, unified set of mosaics 
of the UDF, providing the deepest near-IR blank-field 
view obtained of the universe to date, approaching lim- 
iting magnitudes AB ~ 30 in the near-IR, and yielding 
a legacy dataset on this field that is of lasting scientific 
value. 

2. SURVEY OVERVIEW AND OBSERVATIONAL DESIGN 

This paper presents the overview of the 2012 Hubble 
Ultra Deep Field campaign (UDF12, EST Program ID 
12498, PL: R. Ellis), a large 128-orbit Cycle 19 EST pro- 
gram aimed at extending the previous 2009 WFC3/IR 
observations of the UDF (UDF09, EST Program ID 
11563, PL: G. Illingworth) . The observations were all 
obtained between 4 August 2012 and 16 September 2012, 
and are summarized in Table [TJ The observational ap- 
proach as proposed by the current project is to com- 
bine these and other WFC3/IR imaging observations of 
the UDF area into a single homogeneous dataset, includ- 
ing additional filter wavelength coverage, to provide the 
deepest near-IR observations of the sky currently achiev- 
able, as summarized in Table j4j 

In this project, we aim to study the role of galaxies 
in reionizing the universe, by extending robust searches 
for Lyman-break galaxies to z ~ 9 and beyond, obtaining 

11 http://udfl2.arizona.edu/ 

12 http: / /archive. stsci.edu/prepds/hudfl2/ 



Table 2 

Full-Depth Combined UDF WFC3/IR mosaics (orbits) 



Filter 


UDF09 a 


UDF12 a 


Other a ' b 


Final a 


F105W 


24 


72 


4 


10 


F125W 


34 




5 


39 


F140W 




30 




30 


F160W 


53 


26 


5 


84 



a Number of HST orbits obtained. 

b Other programs include 12099 (PI: A. Riess) and CANDELS 
12060,12061,12062 (PI: S. Faber, H. Ferguson). 



more accurate faint-end luminosity functions at z ~ 7 and 
z ~ 8, and determining more accurate ultraviolet spectral 
energy distributions to constrain stellar populations and 
ionizing photon output. In order to achieve these goals, 
this survey builds upon the previous WFC3/IR invest- 
ment that had already been obtained in the UDF, in a 
number of different ways: 

• Quadruple the exposure time in the F105W filter, 
adding 72 new orbits to the data that had previ- 
ously been obtained in the UDF09 survey (program 
ID 11563, PL: G. Illingworth), in order to provide 
deeper short-wavelength constraints on the z ~ 8 
sources selected using F105W — F125W color crite- 
ria, probe to fainter luminosities, and yield a more 
robust determination of the star formation density 
at 8 -10. 

• Add completely new wavelength information with 
the F140W filter, obtaining 30 orbits of deep in- 
tegration to match the depths in the F125W and 
F160W filters. This provides improved measure- 
ments of the UV slopes of z ~ 7 — 8 sources, ad- 
ditional independent detections of the continuum 
longward of the Lyman break for sources at z > 
9—10, as well as probing to z ~ 12 for sources 
whose Lyman emission may be redshifted out of the 
F140W filter and are detectable only in F160W. 

• Increasing the exposure time in F160W by an ad- 
ditional 26 orbits from the UDF09 program, pro- 
viding more robust red measurements for galaxies 
at z ~ 8 — 10 and further improving the constraints 
on their UV slopes, as well as further securing any 
potential detections up to z ~ 12. 

In this section we provide the details of the observational 
design of the survey, including the filter selection, expo- 
sure times, EST observing considerations, and dithering 
strategies that were employed in obtaining the observa- 
tions. 
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2.1. Filter Wavelength Coverage and Depth 

The science drivers, as summarized above and de- 
scribed here in further detail, drove the filter choice and 
exposure times at a high level to achieve 5a limiting 
depths of F105W = 30.0, F140W = 29.5, and F160W = 
29.5 (AB magnitudes, as measured in ff.'A diameter aper- 
tures), in conjunction with F125W = 29.5 from reprocess- 
ing the previously obtained observations. 

The increase in F105W depth was necessitated by the 
need to improve the robustness of the z ~ 8 dropout se- 
lection as well as probing to fainter luminosities. The 
lack of dynamic range offered in earlier data between the 
F105W — F125W color selections, corresponding only to 
a 2a limit of F105W - F125W > 1.0 at the F125W de- 
tection limits, was inadequate to exclude potential low- 
redshift Balmer-break interlopers. Increasing this dis- 
criminating power to a 2a limit of F105W — F125W > 
1.5, at a ha detection threshold of F125W = 29.5 ob- 
tained from our own reduction of the previous UDF09 
data (program ID 11563, PL: G. Illingworth) , provides 
much stronger selection against low-rcdshift interlopers, 
hence cleaner samples of z > 8 galaxies for studies of 
the luminosity function evolution at these redshifts. We 
achieved this increase with the additional 72 orbits that 
we obtained, thereby reaching our target 5a detection 
limit of F105W = 30.0 (AB magnitudes, 0'.'4 diameter 
aperture) . 

The depth increase for F160W, aiming to match the 
depth of our F125W reduction, was motivated primarily 
by the need to probe further down the luminosity func- 
tion at redshift z ~ 7 and, for the first time, at z ~ 8. In 
particular, establishing whether or not the galaxy popu- 
lation at z ~ 8 can reionize the universe requires measur- 
ing the faint end slope of the luminosity function down to 
jMi500 ~ —17.5, thereby necessitating 5a detection limits 
of ~ 29.5 (AB). This was achieved through increasing the 
existing F160W exposure time by an additional 50% to 
reach a 5a detection limit of F160W = 29.5, from combin- 
ing our new 26 orbits with the previously existing F160W 
data on this field. 

Finally, the addition of the new F140W filter was mo- 
tivated by several science goals. First among these was 
to improve the reliability of any detections of possible 
sources at z ~ 9 — 10 (whose Lyman break has moved out 
of the F125W filter) by using two filters to provide detec- 
tions longward of the break. Since high-redshift galax- 
ies are essentially flat-spectrum sources in /„, this then 
drives the depth of F140W to match that of F160W, i.e., 
a 5a detection threshold of F140W = 29.5 (AB), which 
we achieved using a total exposure time of 30 orbits in 
this filter. This filter also provides additional wavelength 
discriminating power for sources that might be at even 
higher redshift, since Lyman-a would move out of its 
red ward edge at z ~ 12. 

The additional depth in the F105W and F160W fil- 
ters, along with the new wavelength information pro- 
vided by the F140W filter to a matching depth, also 
provide a more accurate measurement of the UV slope 
parameter (3; a robust measure of the average value of 
langle fir angle requires 8a detections to achieve reliable 
constraints. For example, as demonstrated in Dunlop 
et al. (2012), shallower 4a detections can lead to col- 
ors which may be uncertain by up to ~ 0.35 magnitudes, 



Table 3 

WFC3/IR and ACS/WFC Zeropoints 3 - 



Instrument/Camera 


Filter 


Zcropoint (ABmag) 


WFC3/IR 


F105W 


26.269 


WFC3/IR 


F125W b 


26.230 


WFC3/IR 


F140W 


26.452 


WFC3/IR 


F160W 


25.946 


ACS/WFC C 


F814W 


25.947 



a Current information on zeropoints is available at: 

http://www.stsci.edu/hst/acs/analysis/zeropoints 
http: //www. stsci.edu/hst/wtc3/phot_zp_lbn 

b In this paper we present information for F125W, for which we 
did not obtain new observations but instead provide an improved 
reprocessing of existing data, for consistency with the other obser- 
vations that we obtained. 

c Our ACS/WFC data are obtained as parallel exposures on the 
UDF-par2 field, and do not overlap our WFC3/IR observations 
obtained on the main UDF field. 

translating to la errors in j3 of A/3± 1.5. Doubling the 
S/N on the measurements of ft, obtained by using all the 
new data from our program, can therefore estalish much 
more accurate values of langlefi) for galaxies as faint as 
M(/y ~ — 18 at z~7. In addition, the increased depth, 
along with the new F140W information, can extend the 
measured values of (3 for galaxies up to at least z ~ 8. 

In addition to the prime WFC3/IR observations on 
the UDF, we also obtained parallel observations with the 
ACS/WFC camera. We required an orient designed to 
place all these parallel exposures onto the existing par- 
allel field 2 of the UDF (labelled as UDF-PAR2, also 
known previously as UDF-NICP34), using a spacecraft 
orienation of 264° , identical to that which had been used 
for this field in the UDF09 program, to ensure maximal 
overlap with the ACS data that had been obtained as 
part of that program. As discussed further by McLure 
et al. (2011), the existing optical ACS data in that field 
was insufficient to properly exploit the new WFC3/IR 
imaging for the selection of galaxies at z>6.5, and de- 
voting this time to accumulating deep F814W imaging 
for all the parallel exposures was chosen as the most effi- 
cient way of improving this situation for maximum legacy 
value. 

We list in Table [3] the values of the zeropoints cor- 
responding to the four WFC3/IR bandpasses that we 
used. We also indicate locations where more updated 
information may be available, if necessary. These zero- 
points have an accuracy of at least ~ 1 — 2%; remaining 
uncertainties may be related to time-dependent changes 
in the filter or instrument properties, or improved knowl- 
edge of the standard stars that are used in determining 
the calibrations. 

2.2. HST Observations and Dither Patterns 

Due to HST scheduling constraints on our UDF12 pro- 
gram, we divided the 128 orbits into a total of 64 visits, 
where each visit consisted of two orbits, and each or- 
bit consisted of two prime WFC3/IR exposures, accom- 
panied by two parallel exposures using the ACS/WFC 
camera. We adopted similar observing strategies to pre- 
vious programs in order to provide uniformity within the 
final combined datasets. In particular, our dither pat- 
terns followed a strategy consistent with that of the pre- 
vious UDF09 observations (Oesch et al. 2010; Bouwens 
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et al. 2011), which in turn had been based on the orig- 
inal UDF dither strategy (Beckwith et al. 2006). Each 
of the 64 visits (being 4 exposures each) consisted of a 4- 
point dither pattern, shifted around the sub-pixel phase 
space to provide the best possible PSF sub-sampling on 
both the WFC3/IR detector with its relatively large pixel 
scale (tf'128 pixel" 1 ), as well as on the ACS/WFC CCDs 
(0'.'05 pixel" 1 ). In addition, each of these 64 sets of 4- 
point dithers was in turn offset onto a larger-scale grid, 
with offsets of ~ ±3arcsec to cover the ACS/WFC chip 
gap and mitigate large defects and persistence on the 
WFC3/IR detector, as well as providing additional sub- 
pixel phase space sampling. 

All our near-infrared UDF12 observations were ob- 
tained with the WFC3/IR detector, using the IR-FIX 
aperture which samples the full imaging field of view 
(1014x1014 pixels, covering a region ~ 130" across, with 
a plate scale of 0'.'128 pixel" 1 at its central reference 
pixel). All exposures were obtained in MULTIACCUM 
mode using the SPARS100 read-out sequence, enabling 
the pixels to be sampled nondestructively every 100 s. 
These SPARS read samples were repeated for each se- 
quence by specifying the NSAMP parameter, which was 
set to either 15 or 16 before reading out the array, cor- 
responding to either 1300 or 1400 seconds, depending on 
the scheduling constraints for each particular orbit. Each 
of the MULTIACCUM sequences also included two short 
reads at the start, separated by 2.9 s, which served to pro- 
vide a measure of the bias structure across the array at 
the start of each exposure. 

The parallel optical UDF12 observations were all ob- 
tained using the ACS/WFC camera, which comprises 
two CCDs with a usable area of 4096 x 2048 pixels, cov- 
ering ~ 200" in extent with a plate scale of 07 05 pixel" 1 
at the central reference pixel. The two detectors are lo- 
cated adjacent to one another with a small physical gap 
between them of ~ 2'.'5. We used the WFC aperture for 
all the ACS exposures. The exposure times ranged be- 
tween 1200— 1300 seconds, depending upon orbital visi- 
bility constraints. 

3. DATA CALIBRATION AND PROCESSING 

Our final combined UDF09+UDF12 images mo- 
saics have been processed with a version of the 
"MosaicDrizzle" image combination pipeline, specially 
modified for the UDF12 program (see Koekemoer et al. 
2002, 2011 for a more general description). This per- 
forms astrometric alignment and registration, cosmic ray 
rejection, and final combina tion of the exposures u sing 
the MultiDrizzle software (Koeke moer et al7[j2002[ ), as 
well as the Drizzle software ( [Fruchter fe Hook|2002[ ) . In 
this section we provide descriptions ot the input datasets, 
as well as the processing that was carried out within 
"MosaicDrizzle", along with the resulting characteris- 
tics of the mosaics that were produced. We also pro- 
cessed all the UDF09 WFC3/IR data on the main UDF 
in a similar way, as well as all the other overlapping 
WFC3/IR data that were previously described. 

3.1. Initial WFC3/IR Standard Calibration 

We initially processed all our raw WFC3/IR images 
through standard calibration using the Pyraf/STSDAS 



task calwf £p[| in order to obtain a first-pass set of cal- 
ibrated images and carry out initial data quality vali- 
dation. This task populates the bad pixel arrays using 
known bad pixel tables, and subtracts the bias for each 
read using the reference pixels around the border of the 
detector. It then carries out a subtraction of the zeroth 
read in order to remove the bias structure across the de- 
tector, followed by a subtraction of the dark current refer- 
ence files for the SPARS100 read-out sequences. This was 
followed by the non-linearity correction and photometric 
keyword calculation, using the current filter throughput 
tables and detector quantum efficiency curves. 

While the initial calibration was carried out using the 
standard pipeline dark reference files, we found that we 
could improve the signal-to-noise in the final mosaics by 
constructing a custom dark frame from the full set of 
dark calibration files that have been obtained on-orbit 
for the same readout mode and exposure times that we 
were using. Therefore we constructed such a dark frame 
and used it to recalibrate all the exposures, including 
our own as well as those from all the previous WFC3/IR 
observations on this field. 

After having removed basic instrumental effects from 
each read, the exposures were then passed through the 
up-the-ramp slope fitting and cosmic ray rejection steps 
in calwf 3. For each pixel, this step performs a linear 
fit to the accumulating counts that are sampled during 
each MULTIACCUM read, while rejecting outliers from 
the fit as being due to cosmic rays. A final count-rate 
value was then computed for each pixel using only the 
unflagged reads, and was stored as the count-rate in the 
final calibrated exposure, while the uncertainty in the 
slope of counts versus time was stored in the error exten- 
sion of the image. 

3.2. Additional Processing for WFC3/IR Exposures 
with Time- Dependent Variable Sky Background 

A fraction of the UDF12 exposures were affected by sky 
background emission that varied significantly as a func- 
tion of time during the course of the read-out sequences. 
This was most noticeable in the F105W filter, due to 
the wavelength-dependent nature of the background sky 
emission. The resulting time-dependent sky background 
variation introduces a non-linear component into the 
counts that are measured at each read during the MUL- 
TIACCUM sequence. A consequence of this is that the 
standard calwf3 up-the-ramp cosmic ray rejection and 
count-rate slope-fitting algorithm breaks down for these 
exposures, since that algorithm is designed for count- 
rates that are constant in time. Therefore, the resulting 
count-rate images after the completion of calwf3 can 
have significantly non-Gaussian noise properties, which 
can adversely impact the signal-to-noise of the final mo- 
saics, by up to several tenths of a magnitude or more, if 
they are included as-is. 

To mitigate the impact of this issue on the final mosaics 
created for the UDF12 project, we developed a special- 
ized set of routines to carry out intermediate processing 
of the partially calibrated MULTIACCUM sequences, af- 
ter the initial bias and dark current subtraction had been 
performed by calwf 3, but before running the count-rate 

13 Further documentation fo r all the PyRAF/ STSDAS data re- 
duction software is provided at http://stsdas.stsci.edu/ 
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Figure 1. Demonstration of the impact of time-dependent variable sky background on the WFC3/IR images, as well as the correction 
for this effect. All the images are shown to the same greyscale stretch. (Left) initial image after original calibration using the default 
calwf 3, showing strongly non-Gaussian noise properties; (Middle) resulting image once intermediate correction has been applied for the 
time-dependent variable sky and subsequent recalibration, which we have developed for this project; (Right) similar exposure that was not 
affected by time-dependent variable sky; the noise properties of the corrected images (e.g., middle) are now similar to normal images such 
as those on the right, and can be safely combined in the mosaic without adversely impacting its final signal-to-noisc. 



slope-fitting. Specifically, these routines carried out a fit 
to any non-linear components that were present in the 
measured counts at each read, defined in such a way that 
the total integrated non-linear component was zero along 
the time sequence. Thus the non-linear component could 
be subtracted from these reads without any net impact 
on the total accumulated count-rates. 

The corrected exposures were subsequently passed 
through the remaining steps in calwf 3, particularly the 
up-the-ramp cosmic ray rejection and slope-fitting rou- 
tines, which now performed successfully since the incom- 
ing counts were all linear in time. It was verified that 
in all cases, the resulting images had pixel distributions 
that were now perfectly Gaussian, reflecting the mean 
sky background as well as the correct corresponding noise 
properties. Furthermore, the photometry of sources ex- 
tracted from these images was identical to that obtained 
from exposures that had no time-dependent sky back- 
ground. Therefore, after passing all tests on their photo- 
metric and statistical properties, these images were sub- 
sequently included among all the others in the final mo- 
saic combinations, enabling the required depths to be 
achieved. 

3.3. WFC3/IR Persistence, Warm Pixels and Flat 
Field Correction 

After the default calibrations had been applied, we im- 
plemented in our custom pipeline several further correc- 
tions to improve the WFC3/IR data, which are not part 
of the standard calwf 3 pipeline. The first of these con- 
cerns the presence of persistent flux in certain pixels due 
to bright sources having been observed in previous expo- 
sures, which can be a significant issue for the WFC3/IR 
detectorFj In some cases, we were able to make use 
of darks from the WFC3 calibration program that had 
executed just prior to our UDF12 visits, to aid in identi- 
fying and measuring problematic pixels. Pixels with per- 

14 http:/ /www. stsci.edu/hst/wfc3/ins_performance/persistence/ 



sistent flux were then identified in these dark frames if 
they exceeded a count-rate threshold of five sigma above 
the mean, and were flagged in the following science expo- 
sures. For subsequent orbits during a visit, we could then 
determine directly from the preceding exposures which 
pixels may contain sufficient flux to cause persistence; 
the calibration darks were only needed for the first orbit 
in a visit, when the previous data may be from another 
program and not necessarily accessible. In those cases, 
we also made use of the persistence masks created for 
all exposures, accessible from the aforementioned web- 
site maintained by the WFC3 team, which we verified 
were successful in excluding all pixels that were affected 
by persistence. 

We also identified additional "warm" pixels, using the 
full set of on-orbit dark exposures obtained during the 
UDF12 campaign to identify these pixels if they exceeded 
a threshold of five sigma above the mean, in which case 
they were flagged in the data quality arrays that were 
associated with each image, and were excluded from the 
final image combination. We further assembled median 
stacks of all exposures, flagging any pixels that varied 
significantly compared to the general population, which 
resulted in a small number of additional pixels being 
flagged that were not caught from the dark files. 

Finally, the WFC3/IR detector is subject to IR "blobs" 
that have appeared in the WFC3/IR channel since 
launch and were not present in the ground flats, as 
well as some residual large-scale variation in the over- 
all structure of the flatfield. While these are accounted 
for to some extent in the current calibration files, they 
remained noticeable in the deep combined UDF12 imag- 
ing. We therefore made sure to mask out all regions af- 
fected by the blobs in each exposure, as well as applying 
large-scale, low-level residual corrections to the flatfields 
as needed. We verified that the resulting images were 
flat to within ~ 1 — 2% of the mean sky level. 
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3.4. ACS/WFC Calibration 

Each of the raw ACS/WFC exposures in the 
UDF12 parallel field were initially calibrated using the 
Pyraf/STSDAS task calacs. This included bias sub- 
traction, dark current correction, bad pixel masking and 
flatfielding. In addition, a number of other corrections 
need to be applied to ACS data, given the length of time 
that the instrument has been on orbit, as well as account- 
ing for electronic effects in the new CCD Electronix Box 
Replacement (CEB-R) that was installed during Servic- 
ing Mission 4 (SM4). The firs t of these involves t he cor- 
rection for bias striping noise ( Grogin et al.pOlO ) which 
is introduced by the electronics and manifests itself as a 
bias amplitude variation from one row to the next. 

In addition, calacs corrects for the impact of Charge 
Transfer Efficiency (CTE) degradation, whereby charge 
traps present in the pixels can capture some of the elec- 
trons during readout, leading to a loss of flux in the orig- 
inal pixel. This manifests itself as deferred-charge trails 
along the columns behind bright pixels in each exposure, 
while also producing a net astromctric shift up along in 
the column for bright sources. The effect becomes in- 
creasingly severe for pixels furthest from the amplifiers, 
which for these detectors are the pixels near the chip 
gaps. A key point about the CTE correction algorithm 
( Anderson fc Bedin|2010[ ) . is that it is effectively a decon- 
voiution, by virtue of the fact that it restores the charge 
profiles of pixels along a particular column to their origi- 
nal shape, which is sharper and more concentrated than 
the observed profiles which have been smeared by the 
deferred charge trails. As such, the pixel-to-pixel noise 
in the final reconstructed image is also somewhat higher 
than in the original exposure. Tests to date have shown 
that this algorithm correctly reproduces the expected 
noise that would be present in the images if no CTE 
degradation had been present, and that it restores both 
the photometry and the astrometric accuracy to levels 
that are comparable to images without CTE degrada- 
tion. 

Finally, our UDF12 parallel ACS pipelines implement 
a routine to correct for additional bias-related offsets be- 
tween the ACS/WFC3 detector amplifier quadrants, that 
are not fully corrected during standard calibration. This 
routine fits for the differences between quadrants, using 
an iterative clipping procedure to eliminate signal from 
astronomical sources and preserve only the background 
flux, which then removes the residual amplifier quadrant 
differences and places all four quadrants on a uniform 
background level. 

3.5. Relative Astrometry and Distortion 

Once the individual prime WFC3/IR UDF09+UDF12 
exposures and parallel ACS/WFC exposures had all been 
passed through the initial calibrations at the detector 
level, they were subsequently passed through the rest of 
our astrometric and mosaicing processing pipelines. The 
first stage of this pipeline processes all the exposures in 
each visit, for the different cameras, and addresses the 
relative shifts between exposures in each given single- 
orbit visit. The accuracy of the astrometric information 
in the image headers depends on the pointing accuracy of 
the spacecraft, as well as the calibration of the geometric 
distortion models for the detectors. These are described 



in further detail for ACS/WFC in ( |Anderson|[2007| . For 
WFC3/IR we made use of distortion models published 
by (Kozhurina-Platais et al. 2012). 

3.6. Cross- Correlation Shift Determination 

To further solve for and remove the residual uncer- 
tainties in the spacecraft dither offsets between all the 
exposures in each orbit, we applied cross-correlation pro- 
cedures to all the exposures. The cross-correlation proce- 
dure first passes all the exposures through a partial run of 
MultiDrizzle, up to the point where single-drizzled im- 
ages are produced for all the individual exposures, which 
are all aligned on the same pixel grid so that astronomical 
sources should be at the same pixel locations if no resid- 
ual shifts were present. These images were then masked, 
retaining only regions around objects that contain suffi- 
cient signal, and then passed through cross-correlation. 
The resulting cross-correlation peak was fit for each ex- 
posure using a two-dimensional fitting routine to deter- 
mine its location and associated uncertainty, which was 
then directly translated into shifts. This procedure was 
iterated a few times for each exposure, in order to ensure 
convergence of the cross-correlation, with final uncertain- 
ties typically less than a few hundredths of a pixel 

3.7. Final Cosmic Ray and Bad Pixel Rejection 

After having corrected the relative shifts for all the 
exposures, final cosmic ray and bad pixel rejection was 
carried out for all the exposures of a given filter, for each 
camera, by carrying out another run of MultiDrizzle, 
this time with the improved relative shifts. The cosmic 
rays were identified in the driz_cr step of MultiDrizzle 
using a process that first created a series of separately 
drizzled images, one for each input exposure, which were 
subsequently used to create a median image using the 
"minmed" algorithm in MultiDrizzle, which enables 
the minimum to be used instead of the median in cases 
where valid pixels from only two or three exposures are 
present, if one or more of those are impacted by a cosmic 
ray. 

The clean median image was then transformed back 
to the distorted detector frame of each input exposure 
to carry out cosmic ray rejection using the following ap- 
proach. The input counts in a given pixel in the original 
exposure, J e xp, were compared with the counts from the 
median image, / mo d, for the same pixel, together with 
the derivative of the median image, A mec i, defined as the 
steepest gradient from that pixel to its surrounding pix- 
els (with all these quantities being in units of electrons). 
A pixel was flagged as a cosmic ray if it exceeded a cer- 
tain threshold depending on the background sky r.m.s. 
as well as the local gradient around that pixel. The inclu- 
sion of the gradient term A mc d effectively "softens" the 
cosmic ray rejection in regions of relatively steep gradi- 
ents such as bright cores of objects, where the pixel-to- 
pixel variation can exceed simple Poissonian statistics. 
For the data processing in our UDF12 pipeline, this re- 
jection was performed over two iterations, with the first 
pass going through all the pixels in the image and using 
S = 1.2 and SNR = 3.5, followed by a second pass in a 
1-pixel wide region around each of the pixels flagged in 
the first pass, but using more stringent criteria of S = 0.7 
and SNR — 3.0. This ensured that fainter pixels around 
cosmic rays are also flagged. 
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3.8. Absolute Astrometry 

We obtained absolute astrometry for all the visits rel- 
ative to one another by using the original UDF catalogs 
(Beckwith et al. 2006), as well as our own catalog that 
was generated on the UDF-PAR2 field from the exist- 
ing mosaics (Bouwens et al. 2011). All the sources de- 
tected in the exposures for each orbit were matched to 
the sources in the relevant portion of the absolute astro- 
metric catalogs, using a number of iterative steps. The 
first iteration uses a relatively large tolerance (up to a few 
arcseconds) and only the brightest ~ 20 — 30 sources in 
each image, in order to determine the dominant terms in 
the shifts for right ascension and declination. Once these 
had been accounted for, several additional iterations were 
carried out using the full catalog of sources in each image, 
using progressively tighter matching tolerances down to 
071 and solving for the residual remaining shifts as well as 
the rotation errors due to the uncertainties in guidestar 
position. For all visits, ~ 300 — 400 sources were typically 
matched at the faintest levels and tightest tolerances be- 
tween the HST MultiDrizzle-combined images and the 
reference catalogs. 

The results of this procedure are shown in Figures [2] 
andjij for the ACS/WFC and WFC3/IR images. Fig- 
ure Pfshows the distribution in right ascension and decli- 
nation offsets between the reference catalog positions of 
the sources and those measured on the new data, after 
having solved for the astrometry as described above. Fig- 
ure [3] is based on the same data as in Figure [2] but this 
time showing the residuals as a function of position across 
the field, where each vector indicates the mean residual 
in a grid of cells, each 40" on a side, where the size is 
chosen for display purposes to ensure a sufficient num- 
ber of objects per cell while also providing a sufficient 
number of cells across the mosaic to show the general 
structure of the residuals. The residuals are generally 
below < 5 mas for the images, which is comparable to 
the typical on-orbit jitter of HST and sufficient to miti- 
gate all registration uncertainties. 

4. FINAL DATA PRODUCTS 

4.1. Mosaicdrizzle Combination 

We combined all the individual exposures, for each fil- 
ter, into a final mosaics using "inverse variance weight- 
ing", whereby a weight map is created for each expo- 
sure, containing all the uncertainties that are intrinsic 
to a particular exposure (such as dark current and read- 
noise, and background sky noise, modified by the gain of 
the detectors as well as the flat field) . Note that this ex- 
cludes the additional Poisson terms from sources in the 
image, which can be added separately after the fact, if 
needed. 

As an iterative step in the combination, we also re- 
moved remaining low-level background emission in each 
exposure from a masked version of the exposure, which 
was constructed from the full-depth image obtained by 
combining all the individual filters into a single image. 
This image was subsequently smoothed using both a 
small-scale 2-pixel Gaussian, as well as a larger-scale 10- 
pixel Gaussian, which enabled us to construct a mask of 
all the sources on both large scales (to exclude faint outer 
wings) as well as small scales (to exclude small remain- 
ing sources, which would otherwise impact the overall 



sky level determination). 

The pixel scale for the output mosaics is driven by the 
detector plate scale and pixel size, together with the full 
width at half-maximum intensity (FWHM) of the PSF 
produced by the telescope optics. At the wavelengths of 
the WFC3/IR F105W to F160W observations, the HST 
PSF has a FWHM - 0712 - Of! 18, which is subsequently 
convolved by the 07128 WFC3/IR detector pixel scale. 
Hence, the best PSF that could be recovered (without de- 
convolution), even in the ideal scenario of combining im- 
ages using interlacing, which would minimize additional 
convolutions, still has a FWHM ~ 0717 - 0719 in the final 
images. We choose an output pixel scale of 0706 pixel -1 
for the final WFC3/IR mosaics, providing adequate sam- 
pling of the PSF. We also chose a pixfrac parameter of 
0.8, which is small enough to provide some reduction in 
the overall convolution as input pixels are mapped to the 
output plane, while at the same time not being too small, 
so that the overall pixel-to-pixel variation in the output 
weight map is not adversely affected. See Koekemoer et 
al. (2002; 2011) and Koekemoer et al (2011) for further 
details about this parameter and its impacts on the final 
output images, in the context of deep imaging surveys 
with HST. 

Figure [9] shows the full combined mosaics obtained on 
the UDF main field. Three months after the completion 
of the observations, we release all the calibrated mosaics 
to the public via the STScI archiv^J including the driz- 
zled science mosaics as well as the inverse variance weight 
files that describe the noise associated with each pixel. 
Further updates on the project will be provided at the 
primary UDF12 project websit <^®\ as needed. 

4.2. Photometric Limiting Depth Validation 

We have carried out a series of photometric and limit- 
ing depth tests on the full combined WFC3/IR mosaics, 
aiming to validate the depth achieved in absolute terms 
as well as relative to the previous data on this field. In or- 
der to quantify the limiting depth across the mosaic, we 
first constructed a full-depth, full-filter mosaic using all 
the WFC3/IR observations on the UDF, in all four filters 
(F105W, F125W, F140W, F160W). This broad-band im- 
age provided extremely deep sensitivity for masking out 
the extended faint wings of sources, which is necessary in 
order to obtain genuinely source-free regions in order for 
an accurate background r.m.s. estimate to be obtained. 
In order to improve the signal-to-noise with which faint 
objects were masked, we created smoothed versions of 
this full-depth image by convolving it with a small-scale 
Gaussian (2 pixels FWHM, aimed at detecting and mask- 
ing all the faint sources), as well as a larger-scale Gaus- 
sian (10 pixels FWHM), which successfully masked all 
extended emission around larger sources. 

In Figure[4]we show the full-depth image obtained from 
all four filters, which was used to create the object mask 
that we subsequently applied. The object mask excludes 
about 45% of all the pixels in the mosaic; the remain- 
ing pixels were then considered to represent the pure 
sky background (along with potentially exceedingly faint 
sources that are not included in the object mask). The 
statistics of these pixels were analyzed using several dif- 

15 http:/ /archive. stsci.edu/prepds/hudfl2/ 

16 http://udfl2.arizona.edu/ 
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Figure 2. Residual uncertainties in astrometry for a total of ~ 1400 objects detected in the main UDF12 field, that were cross-correlated 
with sources from the original UDF catalog by Beckwith et al. (2006), which was used for the astrometric registration for this project. 
Each point on the plot represents the remaining positional offsets between the sources, after having corrected for all the global < 0'.'005, 
after accounting for the positional uncertainty of each of the sources. As a result, the astrometry is sufficiently accurate to allow all the 
exposures in the dataset to be reliably combined, with no significant residuals remaining. 
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Figure 3. Residual uncertainties in astrometry represented as vectors on the UDF field, for each of the individual UDF12 sources that 
were matched to the original UDF catalog by Beckwith ct al. (2006). This indicates that there are no regions of net residual astrometric 
shift remaining across the image, thereby further confirming that the accuracy of the overall alignment is robust to the level of < 0'.'005, 
and sufficient to enable reliable combination of the full set of input exposures. 
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Table 4 

Final limiting 5cr sensitivities' 1 



Filter 


Negative Gaussian 


Average r.m.s. 


Global r.m.s. 


F105W 


29.98 


29.97 


29.97 


F125W 


29.55 


29.53 


29.52 


F140W 


29.51 


29.49 


29.48 


F160W 


29.46 


29.45 


29.45 



a All measured depths are for our final, full combined mosaics, 
representing the 5cr limiting depth in apertures of diameter 0'.'4, 
Of/ 44, C//47 and Of.' 50 respectively for F105W, F125W, F140W and 
F160W. 



ferent tests, in order to determine the r.m.s. values on 
small scales as well as determining the global uniformity 
of r.m.s depth across the mosaic. These pixel-to-pixel 
r.m.s. values, after accounting for correlated noise, then 
provide a direct estimate of the limiting sensitivity of the 
mosaic. 

The first test involved dividing the mosaic into a regu- 
lar grid of cells, in order to determine the relative degree 
of depth variation across the mosaic, as well as the degree 
of flatness in the residual sky level. This test reveals the 
degree of uniformity in detection sensitivity between dif- 
ferent cells, as well as the impact of any remaining large- 
scale residuals in the background sky, which would serve 
to broaden the global measured r.m.s. as compared with 
the average of all the individual r.m.s. measurements 
obtained in the different cells. The results from these 
tests are shown in Figures [HI— [8j where we translate the 
pixel-to-pixel r.m.s. into 5crlimiting magnitudes in aper- 
tures of diameter C//4, f/.'44, Of.' 47 and Of! 50 respectively 
for F105W, F125W, F140W and F160W (corresponding 
to 70% of the total enclosed flux), after having accounted 
also for the presence of correlated noise by comparing the 
science images with the inverse variance weight maps. 

In addition, the limiting depths are relatively uniform 
across the field, with residual variations <0.03 mag- 
nitude except for occasional areas of somewhat lower 
weight that correspond to known large defects on the de- 
tector. Finally, the global measured r.m.s. values from 
these tests agree very well with the average r.m.s. from 
the individual cells (<0.01 magnitudes), indicating that 
the background residual sky structure is globally flat with 
no significant impact on the overall r.m.s. These results 
are presented in Table HI 

We also performed blank-aperture tests, calculating 
the 5cr depths based on a total of 15,000 apertures placed 
within blank sky regions across the mosaic. These tests 
include fitting only the negative half of the pixel dis- 
tribution, in order to exclude low-level positive sources, 
and used the same sized apertures as the previous tests. 
These results are also presented in Table [4] The re- 
sults agree very well, indicating that the depths achieved 
match the expected sensitivities for each of the filters, for 
these full-depth mosaics, demonstrating that we achieve 
our proposed limiting depths of AB ~ 30 for F105W, and 
AB~29.5 for F125W, F140W and F160W. 

5. SUMMARY 

We have described the 2012 Hubble Ultra Deep Field 
campaign (UDF12), a large 128-orbit Cycle 19 HST pro- 
gram aimed at extending previous WFC3/IR observa- 
tions of the UDF by quadrupling the exposure time in 



the F105W filter, adding a completely new F140W fil- 
ter, and extending the F160W exposure time by 50%. 
The project is aimed at determining the role played by 
galaxies in reionizing the universe, and includes obtain- 
ing a robust determination of the star formation density 
at z>8, improving measurements of the ultraviolet con- 
tinuum slope at z r~j 7 — 8, facilitating the construction 
of new samples of z ~ 9 — 10 candidates, and enabling 
the detection of sources up to z~12. For this project 
we committed to combining these and other WFC3/IR 
imaging observations of the UDF area into a single ho- 
mogeneous dataset, to provide the deepest near- infrared 
observations of the sky currently achievable. We have 
described the observational aspects of the survey as mo- 
tivated by its scientific goals, and have presented a de- 
tailed description of the data reduction procedures and 
products from the survey. We release the full combined 
mosaics, comprising a single, unified set of mosaics of 
the UDF, providing the deepest near-infrared blank-field 
view obtained of the universe to date, reaching magni- 
tudes as deep as AB ~ 30 in the near-infrared, and yield- 
ing a legacy dataset on this field of lasting scientific value 
to the community. 
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Figure 4. Full-depth UDF image, created by combining all the F105W, F125W, F140W and F160W exposures, that was used to create 
the object mask subsequently used for the blank sky statistical measurements for the limiting depth calculations. 
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Figure 5. Measured 5cr limiting sensitivity across the full-depth F105W UDF mosaic. See text for further details. 
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Figure 7. Measured 5<r limiting sensitivity across the full-depth F140W UDF mosaic. See text for further details. 
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Figure 9. Image showing the full mosaics from our new WFC3/IR UDF12 data combined with the previous UDF09 and other data on 
this field, including all the WFC3/IR filters (F105W, F125W, F140W, F160W). 



